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Abstract 



We investigate a microscopical structure in a chain of cars waiting at a red signal on signal-controlled crossroads. Pre- 
sented is an one-dimensional space-continuous thermodynamical model leading to an excellent agreement with the data 
. measured. Moreover, we demonstrate that an inter-vehicle spacing distribution disclosed in relevant traffic data agrees 

' with the thermal-balance distribution of particles in the thermodynamical traffic gas (discussed in UJ) with a high inverse 

temperature (corresponding to a strong traffic congestion). Therefore, as we affirm, such a system of stationary cars can 
C"| ■ be understood as a specific state of the traffic sample operating inside a congested traffic stream. 

Key words: vehicular traffic, parking problem, particle gas, spacing distribution, thermal equilibrium, Random Matrix 
Theory 
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C\ . 1 Introduction 

, Investigation of various transport systems is recently one of the prominent subjects of physics. Intention of relevant re- 

t— ( ' searches is to describe such systems (or phenomena) quantitatively, create their appropriate models (theoretical or nume- 

C**"* . rical), and finally obtain the exact or numerical outputs comparable to real situations. Higher aspiration of those scientific 

' researches can be found in finding a certain connection among the different phenomena and revealing a possible univer- 

, sality. 

•rH . 

Currently, one of the strongly accented fields is an investigation of queuing systems. Under the terms of that field it 
C$ ' has been discussed many various topics, for example, wide-ranging spectrum of vehicular traffic problems [2 |, pedestrian 
dynamics Q, escape panic flU, longitudinal parking of cars on a street JB |, [6|, parallel parking J7J, JD, or public transport 
in some Latin America countries (9), iflOl . All those subjects are in close connection with the Random Matrix Theory, 
theory of chaos, or theory of particle gases (see the references cited above). The main goal of this paper is to extend the 
set of queuing systems mentioned above by a stationary ensemble of cars waiting at a red signal on signal-controlled 
crossroads (see also ITD ). 

Moreover, we are aiming to create an one-dimensional model of point-like vehicles producing the same inter-vehicle 
gap distributions as those detected among cars standing on signal-controlled intersections. In the second part of this article 
we demonstrate that such a model can be interpreted (on a microscopical level) as a thermodynamical gas of dimensionless 
particles exposed to a thermal bath. This analogy allows, as we assert, to find an exact form of relevant spacing distribution 
which can be consequently compared to the realistic gap-statistics. 



2 Describing the system 

The traffic data analyzed in this work were measured during a few days on a multi-lane intersection located near center 
of Prague. This intersection is a constituent of an extensive network of roads and crossroads inside the internal metropolis 

*e-mail address: milan . krbalek@ f j f i . cvut . cz 
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Fi gure 1 ; Graphic representation of the model. The upper 
subplot depicts the initial state of the numerical scheme de- 
scribed in text whereas the bottom subplot demonstrates the 
final stationary state of the traffic sample, i.e the state when 
the cars are waiting for green signal. We note that the squares 
represent the model particles with the leading car being picked 
out. 



and is therefore strongly saturated during the whole daytime practically. Furthermore, the time interval between two green 
signals (on one crossroad) is very short, which causes that some cars are not able to reach the threshold of following inter- 
section (during one green phase) and have to wait therefore for another green light. This fact finally leads to a substantial 
decrease of average velocity of vehicles moving between crossroads, i.e one can observe the effects usually detected in 
congested traffic regime (see Ref. |[2lD. Measured were bumper-to-bumper distances between subsequent cars ((i + l)th 
and ith ones) waiting at a red signal (in one direction only). Data file contains 5022 digitally gauged events showing the 
mean inter-vehicle gap equal approximately to 149 centimeters. We add that the clearances were measured directly using 
the laser technology. 

More detailed statistical analysis uncovers that a probability density p(r) for distance r between neighboring cars 
shows a similar behavior as that investigated between the eigenvalues of random matrices (see |[T2l ). zeros of Riemann 
zeta function (see lfT3l ). or vehicles moving inside the traffic stream on the freeways (see 0]). Such a behavior (follow 
Fig. 2) demonstrates the presence of repulsive interactions among the elements in question. As well known, a spacing 
distribution of non-interacting elements shows a different distribution, in concrete: Poisson probability density 

p(r) = cxp[— r] (r > 0). 

Since the traffic interaction (in local sense, of course) is usually quantified as power-law repulsion among the successive 
vehicles (see Ref. (T) and lfT4l ) let us suppose that a potential energy of the ensemble investigated reads as 

N 

U(r u r 2 ,...,r N )=Y / r- 1 - (D 

i=l 

Herein we assume that the stationary traffic state analyzed in this paper (i.e. the queue of waiting cars) is determined by 
the preceding process - traffic flow towards the intersection. Evidently, moving in traffic sample the driver is interacting 
with other cars and optimizing his/her motion to reach the threshold of the crossroad as soon as possible and, at the same 
time, avoid a crash with the preceding vehicle. Such a behavior corresponds to the thermodynamic effects governing the 
ensemble into a local thermal equilibrium (see JT] for details). 



3 Modified Metropolis algorithm 

Accepting the above-mentioned assumptions on thermodynamical aspects of the issue we formulate the following one- 
dimensional traffic model based on principles of statistical physics. Consider N + 1 point-like particles (cars) located 
randomly (or equidistantly if advantageous) on a line (or on a circle) so that the mean gap among them is one, i.e. 

N 

J>=iV, (2) 

i=X 

where r, represents the gap between (i + l)th and ith particles. Thus, the ordered positions x\ > x 2 > ■ ■ ■ > a;jv+i 
constitute the initial state for our simulation (see Fig. 1). The particles move along the line (or along the circle) accepting 
the undermentioned rules until the leading car reaches a fixed point (the threshold of new crossroads). In accord to a 
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Figure 2: Inter-vehicle gap statistics p(r). Bars represent 
the probability density for bumper-to-bumper distance among 
the cars waiting at a red signal on intersections (measured 
in Prague). Data were re-scaled so that the mean spacing is 
equal to one. Points display the optimized result of the nu- 
merical scheme (Metropolis algorithm) for N = 100 and 
Pmodd — 1-45. Finally, the curve displays the distribution 
{3) for the fitted value f3f it ss 1.2488 (obtained by the num- 
ber variance test.) 

realistic situation the overtaking cars are not permitted, i.e. the particles can not change their order. Let fi mo dei > denote 
the inverse temperature specifying the measure of chaos inside the ensemble simulated. We assume (3 mo dei to be the only 
significant parameter of the model. The car positions are repeatedly updated (we use 20000 steps in our version) according 
to the following rules: 

1. Calculated is the potential energy Uq (using formula (HJ) for the actual set of locations {x\, X2, ■ ■ ■ , xn+i}- 

2. We pick an index j £ {1, 2, . . . , TV + 1} at random. 

3. We draw a random number 8 equally distributed in the interval (0, 1). 

4. We compute an anticipated position x'j = Xj + S of jth element. Because of singularity in the potential energy ((TJ 
the model particles can not change their order. Therefore we accept x'j only if x'j < Xj-\. 

5. We calculate a value of potential energy U' determined for configuration {xi, X2, ■ ■ ■ , Xj—i,x'j, Xj+i, . . . , xjv+i}- 

6. If U' < Uq the jth particle position take on a new value x'j . If V > Uq then the Boltzmann factor 

w = cxp [-/3 model AU] , 

where AU = U' — Uq, should be compared with a random number r equally distributed in (0, 1). Provided that 
the inequality w > r is fulfilled the jth particle position takes on the new value x'j too. Otherwise, the original 
configuration {xi, X2, ■ ■ ■ , xn+i} remains unchanged. 

The sketched procedure represents a modified Metropolis algorithm originally developed for chemistry purposes (in 
Ref. |15|). This algorithm belongs to the category of Monte Carlo simulations (see Ref. ifTTl ) which are recently used 
for numerical modelling of statistical systems (as demonstrated in Ref. JT3), for example). The elaborated scheme of 
Metropolis ensures a relaxation of ensemble into a thermal-balance state when the energy fluctuates around a constant 
value being independent of initial configuration of particles (see Fig. 3). After reaching the thermal equilibrium (i.e. after 
approximately 5000 updates of configuration (Monte Carlo steps), as visible in Fig. 3) the ensemble lingers in this state 
until the simulation is interrupted. Then, as observed, corresponding probability density for inter-particle gaps depends on 
the inverse temperature (3 m odei only. 

Our aim is to find the optimal value of inverse temperature (3 m odei so that the gap distribution p(r) corresponds to that 
measured among the cars on crossroads. Using a x 2 — method (i.e. minimizing the sum of squares-deviations between two 
distributions in question) one can find that optimal value (3 mo d e i is approximately 1.45. Concretely, for fixed value (3 mo dei 
the distribution p(r) is obtained. Then the x 2 — test between empirical data and p(r) could be evaluated. The optimal 
value of Pmodei is the one for which the corresponding sum of squares-deviations is minimal. To conclude, for value 
Pmodei = 1-45 both processes (traffic and Metropolis procedure) generate practically the same gap distributions (see Fig. 
2). Thus, the introduced procedure could represent a realistic model for behavior of the cars in the vicinity of the chosen 
intersection. 
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Figure 3: Relaxation of the system into the thermal equi- 
librium. Dashed and continuous lines (see the upper-left or 
lower-left corners, respectively) display the energy value {T} 
for N = 100 and /3 mo( j e ; = 1.45 during the run of Metropo- 
lis procedure (having 20000 steps) for random (or equidistant) 
initial locations of elements, respectively. Plotted is the aver- 
age value (calculated for 100 repeated realizations of Metropo- 
lis). Grey curve represents the energy value {T} for one realiza- 
tion of Metropolis (when initial particle positions were chosen 
equidistantly). 



4 Terminal state of thermodynamical traffic gas 

As explored in articles lfl4l . JTj, and ifTHl . the traffic flow can be understood (on a microscopical level) as a thermodynam- 
ical gas of interacting cars exposed to a heat bath of inverse temperature (3. Besides, the latter has an immediate relation 
to the traffic density. If accepting such a approach we describe the traffic ensemble (on the move) as a circular gas of 
point-like particles whose hamiltonian reads as 

i i 

where Vi and r, represent a ith car velocity and gap to the previous car, respectively. Quantity v denotes the desired 
velocity of the ensemble. Then (see the exact calculation in [[Q) the derived probability density pp(r) for a gap r among 
the successive vehicles is 

pp(r) = Aexp^r- 1 -Br] , (3) 
where the constants A and B are calculated via two normalization equations 

r pp(r) dr = 1. 



pp{r) dr = 

According to Ref. Q] the following relations hold true: 

B«/3 + 



\/2/3 + 3-e-v^ 



V&filCi (y4/3 2 + 6/3 - 2/3e-V73j 
Herein /Ci stands for a Mac-Donald's function (modified Bessel's function of the second kind) of the first order. 

Since the situation investigated in this article is without any doubt the result of a preceding traffic flow (see ifTTI ) it 
is meaningful to expect that the clearance distribution among the cars waiting at the red-light-signal will be of the form 
(0). Indeed, as confirmed by an appropriate statistical analysis of the collected data (discussed later) the measured gap 
statistics (see Fig. 2) corresponds to the probability density (O if the inverse temperature (3 of the thermodynamical model 
is 

p fit w 1.2488. (4) 

We denote that this value has been determined by a more sophisticated method presented in next section. In addition, 
a positive comparison between the corresponding gap distributions supports the hypothesis that traffic stream can be 
locally understood as a stochastic gas whose elements are repulsed by the forwardly-directed nearest-neighbor power-law 
potential depending on a reciprocal distance between successive gas-elements. This correspondence, however, does not 
mean that traffic is a thermodynamical system, of course. 
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Figure 4: Results of the number-variance test. Plus signs 
stand for the value A n (L) calculated for the collected traffic 
data, whereas points represent the same quantity for particles 
of Metropolis model (where f3 mode i = 1.45). Dash-dotted 
line visualizes function A n (L) = L representing the number 
variance of independent events. The solid curve displays the 
function j5J calculated for the optimal value ft fit ~ 1.2488 
obtained by the X 2 — fit to the traffic data. 



5 Testing the statistical variance of data 

If trying to find a more robust argumentation for an assertion on statistical similarities between the process investigated 
and the traffic model we can apply some of the techniques originally developed for purposes of the Random Matrix The- 
ory (see the book lfl2l ). Usual way how to quantify the behavior of variances among the statistical data is in applying 
so-called number-variance test. Such a test is defined as follows. 

Consider TV spacings n, r2, . . . , rjv between the successive vehicles (or particles of model) and suppose that the mean 
distance taken over the complete ensemble is re-scaled to one, i.e. 

N 



J2n = N. 



2 = 1 



Dividing the interval [0, TV] into subintervals \{k — 1)L, kL] of a length L and denoting by rik(L) the number of cars in 
the fcth subinterval, the average value n(L) taken over all possible subintervals is 

x VN/L\ 

= W/L} S nk(L) = L ' 

where the integer part \_N/L\ stands for the number of all subintervals included in the interval [0, TV]. Number variance 
A n (L) is then defined as 

and represents (in a traffic instance) the statistical variance in the number of vehicles operating at the same time inside a 
fixed part of the road of a length L. 

As well known from Random Matrix Theory the number variance can be explicitly derived from the relevant spacing 
density pp{r). The significant advantage is remarkable sensitivity of the number variance A n (L) to any change in the 
probability density ppir) - i.e. to any change in the potential U(ri,r2, ■ ■ ■ ,r n ) also. Whereas the number variance of 
independent events (or non-interacting elements) is the identity A„(L) = L, for a thermodynamical traffic gas with non- 
zero inverse temperature j3 there has been numerically calculated (in Ref. Ifl8l ) a different behavior, concretely: a linear 
dependence 

A n (L)^ x L + j (5) 

with a slope 

1 

2.4360 f^ 8 ^ 

and a shift 

P 



Y PS < 1 

A i Aoan <30.8207 .11 — 



5.1926/3 + 2.3929 



> 0. 
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As understandable now, the comparison between the number variance of the collected data and the function (f5]) can be 
then used (together with the comparison of the relevant gap distributions) as a robust fitting procedure which is capable 
of revealing more detailed nuances among the distributions compared. If applied to our topic, such a procedure generates 
the optimal value (0]i for which the exactly determined number variance (0) corresponds to the measured data (see Fig. 
4). Note that both of curves A n (L) are rapidly deflected from the line visualizing the number variance of non-interacting 
particles. It implies the presence of a strong repulsion among the vehicles. However, a small deviation is detected for 
larger L between the traffic data (plus signs in Fig. 4) and Metropolis data (points in the same figure). Such a discrepancy 
can be explained by the simple fact that the respective temperatures (i.e. [3 mo dei and /3/it) differ each from other. 

6 Summary and discussion 

Investigated was the traffic ensemble of vehicles waiting at a red-light-signal on signal-controlled crossroads. We have in- 
troduced the thermal space-continuous time-discrete traffic model of repulsing point-like elements based on the Metropolis 
algorithm. By the suitable choice of the inverse temperature parameter there were obtained the same statistical distribu- 
tions as those produced by the real traffic process. Above that, we show that the investigated state of the realistic traffic 
sample can be predicted with the help of the thermal-equilibrium state for local thermodynamical gas whose point-like 
particles are repulsed by the short-range power-law potential (Q]). As demonstrated above, for the fitted value j3 ~ 1.25 of 
reciprocal temperature the corresponding spacing distributions are practically the same. The correspondence between the 
traffic samples and presented theory is, moreover, supported by the robust test of number variance which reveals 

1 . the thermal feature of the topic - on microscopic scale, 

2. the presence of strong interactions among the cars, 

3. a deep connection between the stationary state of waiting cars and the preceding move of sample towards the 
intersection threshold. 

To conclude, we assert that the configuration of vehicles waiting at a red-light-signal on signal-controlled crossroads is 
a product of local thermodynamics-like processes acting among the cars. All the accessible statistical analyses strongly 
support this fact. Therefore, the observed phenomenon can be understood as a traffic in especial super-congested state. 
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